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A high-resolution neutron spectroscopic technique is used to measure momentum-resolved magnon 
lifetimes in the prototypical two- and three-dimensional antiferromagnets Rb2MnF4 and MnF2, over 
the fuh Brillouin zone and a wide range of temperatures. We rederived theories of the lifetime result- 
ing from magnon-magnon scattering, thereby broadening their applicability beyond asymptotically 
small regions of wavevector and temperature. Corresponding computations, combined with a small 
contribution reflecting collisions with domain boundaries, yield excellent quantitative agreement 
with the data. 
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Dynamics in Heisenberg antiferromagnets have been 
studied very extensively by means of theoretical ap- 
proaches such as hydrodynamics [T], dynamical scaling[5], 
and interacting spin- wave methods [3]. In the dilute 
magnon gas in two and three dimensions [IHl], the dom- 
inant damping processes are identified by diagrammatic 
perturbation theory to be elastic pair-wise collisions of 
magnons (see Fig. la-b); the transition rate is given by 
Fermi's Golden Rule. Analytic formulae for the relax- 
ation as a function of wavevector q and temperature T 
were derived for this mechanism (see Fig. 1). Follow- 
ing the discovery of high-temperature superconductivity 
in doped two-dimensional (2D) antiferromagnets, recent 
theoretical work has focused specifically on the differ- 
ences between the behavior of magnons in two and three 
dimensions E HO]. Surprisingly, in recent experimen- 
tal studies of both 2D [11] and 3D [12] antiferromagnets, 
the measured linewidths disagreed considerably with pre- 
dicted values, and anomalous additional line broaden- 
ing at low temperatures that eluded explanation was ob- 
served. To address these issues, we have reexamined this 
decades-old problem both experimentally and theoreti- 
cally, in both 2D and 3D. 

Rb9MnFA[T5] and MnF2[Tl] are near-ideal realizations 
of 2D and 3D antiferromagnets, respectively. Due to 
their structural and chemical similarity, they are a well- 
suited pair of compounds for comparison of linewidths in 
2D and 3D. In these body-centered tetragonal systems, 
the magnetic interactions are dominated by the antifer- 
romagnetic (AF) exchange coupling J between 5 = 5/2 
Mn^+ spins: for Rb2MnF4, between the 4 nearest neigh- 
bors in the a6-plane, and for MnF2, between the 8 next- 
nearest neighbors (arranged at the corners of a tetragonal 
prism). A uniaxial anisotropy parameter D arising pri- 
marily from dipole-dipole interactions causes the spins 



in both compounds to align along the c-axis. The re- 
sulting dispersion relations in the two systems are very 
similar: the zone-boundary energy and energy gap at 4 
K are 6.6 meV (along (qqO)) and 0.6 meV in Rb2MnF4, 
respectively |T3), and 6.3 meV (along (gOO)) and 1.1 meV 
for MnF2[T5]. The Neel temperature Tn is 38.4 K in 
Rb2MnF4 and 67.6 K in MnF2. 

Wc used the high-resolution NRSE-TAS (neutron res- 
onance spin-echo triple- axis spectroscopy) technique |16j 
to measure magnon linewidths in the a6-plane of these 
antiferromagnets (see Figs. 2-4). Temperatures ranged 
from 3 K up to 0.8 Tn (Rb2MnF4) and 0.6 Tn (MnFa); 
the wavevectors spanned the full magnetic Brillouin zone 
(BZ). The neutron spectrometer and experimental meth- 
ods used are as described in Ref. 12. The mosaicity of 
the samples used was measured by means of gamma-ray 
diffractometry and found to be 0.44' and 0.99' for MnF2 
and Rb2MnF4, respectively. The MnF2 sample was ori- 
ented in the (HKO) scattering plane and the Rb2MnF4 
sample in the (HHL) scattering plane. The magnetic BZ 
used in MnF2 was that centered around (010) and that in 
Rb2MnF4 around (HO) (the latter in tetragonal crystal 
coordinates with a — 4.216A; this reflects K2NiF4-type 
magnetic ordering|13j). The final neutron wavevectors kf 
used ranged from 1.7 to 3.2 A ^ . 

In NRSE-TAS magnon linewidth measurements, accu- 
rate knowledge of the magnon dispersion is required in or- 
der to set the angles of the RF tilt coils correctlv |12[|17) . 
Correspondingly, we measured the dispersion at all tem- 
peratures at which we measured linewidths (except where 
the change with temperature was not experimentally re- 
solvable) . 

The results for both compounds are qualitatively sim- 
ilar to those found for gc-magnons in MnF2 [121 • A broad 
peak at 3 K of maximum value ~6 fj,eV at small q charac- 
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FIG. 1. (Color online) (a),(b) The primary scattering pro- 
cesses involve four magnons. (a) The two-in/two-out process 
dominates at low temperatures, (b) The one-in/three-out (or 
three-in/one-out) process is calculated to be negligible com- 
pared to that of (a), (c) The dispersion surface in 2D. The 
lifetime broadening is determined by the number of accessi- 
ble final states. These scattered states lie on loops on the 
surface, shown here for k = (0.1,0.1) and p = (0.05,0.2). 
The loop is also projected onto the Brillouin zone. Strict con- 
servation of energy and momentum results in the loop being a 
delta function; lifetime broadening relaxes this condition. For 
small (large) k, the loop shrinks (grows), resulting in narrow 
(few scattering states available) and broad (many available) 
linewidths, approximately as ~ fc. In 3D, the scattered 
states are in the surfaces, or "shells", of 3D volumes, with 
TsD ~ k^- (d) For illustrative purposes, we have concretized 
the definition of the four regimes for which TH provide ana- 
lytical solutions in 2D by stipulating that a strict inequality 
represents a factor of 10 difference in magnitude. The cor- 
responding boundaries of the regimes A-D [5j are indicated. 
The grey line indicates the highest temperature at which we 
measured linewidths. The region of applicability of Kopietz' 
analytical result ^ is shown in orange; the upper limit of 
the temperature does not quite overlap with the lowest tem- 
perature at which we measured (3 K). Analogous regimes are 
employed by HKHH in 3D; their area is comparable to that 
shown for TH. 



terizes both data sets. We initially compared the results 
with analytical approximations for the magnon linewidth 
originating from two-in/two-out magnon collisions (see 
Fig. la). The corresponding expressions derived by Har- 
ris et al. (HKHH) 4] and Tyc and Halperin (TH) [5 ap- 



ply to restricted regions of small wavevectors at low tem- 
peratures (see Fig. Id). Their general low-temperature 
expressions, which contain no adjustable parameters and 
involve summations over the BZ, are intended for use at 
long wavelengths. Umklapp processes, which become in- 
creasingly important as q and T increase, are not treated 
properly in the scattering matrix elements. 

We rederived the scattering matrix elements for gen- 
eral q, treating Umklapp processes correctly [TBj, and 
then computed the corresponding linewidths rniag((7, T) 
numerically. It became clear that "on-shell" computa- 
tions do not give correct line broadenings, particularly 
near the zone boundaries, because of the curvature of the 
dispersion relation and the large density of states present 
there. For this reason, we evaluated Tn^agiq, T) self- 
consistently. We incorporated experimentally-measured 
values of the renormalized magnon energies into the cal- 
culation, and hence our results can be applied at higher 
temperatures than those considered by HKHH and TH. 
Through such calculations, we determined that incor- 
rect treatment of Umklapp processes yields inaccurate 
linewidth results, even at small wavevectors and low tem- 
peratures. Independent of how these are treated, our 
numerical results differ considerably in magnitude from 
the analytical approximations, and thus numerical eval- 
uations are indispensable even for the limiting cases of 
small q and low T. 

The computations involve sums over p and q states, 
where k-|-p = r-fs, r = p — q, and s = k -|- q. For p, 
a grid with density of {81 / {2tt))'^ pts ■ a?, with 64 times 
higher density in the central 1/8*'' of the BZ, is used 
in the summation, weighted by the reciprocal space vol- 
ume per point. The grid in q has a uniform density 
of (81/(27r))^pts • . Different grid sizes were checked 
to ensure that convergence occurred. Momentum con- 
servation in the scattering process is achieved by the 
choices r = p — q and s = k -I- q. To approximate the 
delta-function energy-acceptance window, a Lorentzian 
probability distribution p(A) = l/(7r7[l -I- (A/7)^]) is 
used, where A = Lij{k) + a;(p) — a;(r) — u}{s) is the devia- 
tion from the conserved result, and the width 7 is deter- 
mined by the lifetimes of the "intermediate" spin-waves 
involved in the scattering process, viz. 7 = 7p + 7r + 7s 
where the intermediate line broadenings are interpo- 
lated for the relevant wavevectors from those previously 
computed. To speed up the calculation, the Lorentzian 
is truncated at ±47, with its amplitude adjusted cor- 
respondingly to preserve unit area. Self-consistent cal- 
culations are performed by choosing an initial narrow 
linewidth that is uniform in wavevector and then iterat- 
ing. Iterations were performed until the results converged 
to within 0.4% (for MnFz) or 1% (for Rb2MnF4), over 
the full range of wavevector and temperature. 

An additional potential scattering process consists of 
a non-momentum-conserving collision of a magnon with 
a boundary (e.g. a grain boundary, antiferromagnetic 
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FIG. 2. (Color online) Linewidth data over the full BZ in 
(a) Rb2MnF4 and (b) MnF2. The curves show self-consistent 
calculations, with no free parameters, of 4-magnon scattering, 
added to a small, nearly T-independent contribution, peaked 
at low q, representing 2-magnon scattering from crystal blocks 
(with one free parameter, fitted to the 3 K data). 



domain wall, or microcrack in the crystal). A rough 
estimate of the associated magnon lifetime for such a 
process is = L/{2v{q)), where L is the size of the 
grain, domain, or microcrystallite and v{q) is the magnon 
velocity |19j. The corresponding additive contribution 
Fb = hr^^ to the overall linewidth exhibits a peak at 
low q that is similar to that in the 3 K data (see Fig. 
3). We fitted these data to the corresponding expression 
after first subtracting off the calculated Fmag- Using the 
measured magnon dispersion data, we obtained an aver- 
age L of 0.58 ± 0.04 ^im for Rb2MnF4 and 0.50 ± 0.04 
/im for MnF2. The latter is consistent with the antifer- 
romagnetic domain size of 0.63 ± 0.05 /im obtained from 
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FIG. 3. (Color onhne) Linewidths at 3 K for Rb2MnF4 (red) 
and MnF2 (blue), with the calculated 4-magnon scattering 
subtracted off (this is less than 0.4 fieV for Rb2MnF4 and 
0.05 fieV for MnF2, for all q). The dashed curves show fits 
to 2-magnon scattering (see text); the fits shown by the solid 
curves also allow for a constant offset. 



independent high-resolution neutron Larmor diffraction 
measurements performed on the same MnF2 crystal at 5 
K. 

Under the assumption that L remains constant at all 
temperatures, we calculated Fb for higher temperatures 
(again using the measured magnon dispersion data). The 
corresponding temperature dependence is small: less 
than 10% of the 3 K values, over all q, at the highest 
temperatures measured for MnF2, and less than 5% for 
Rb2MnF4. 

The fits to the 3 K data can be improved considerably 
by adding a small, temperature-independent constant (or 
such a constant plus a positive, temperature-independent 
slope) to the above contribution from 2-magnon scatter- 
ing (see Fig. 3). This constant is found to be 2.1 ± 0.4 
HeV for MnF2 and 1.7 ±0.3 /leV for Rb2MnF4. The cor- 
responding values of L are 1.4 ± 0.8 fim for MnF2 and 
1.0 ± 0.2 /im for Rb2MnF4. Such an offset could refiect 
either an additional source of linewidth broadening or 
a systematic experimental error of as-yet undetermined 
origin. 

One mechanism that could cause a line broadening 
of this order of magnitude is scattering from the finite 
width of the crystal boundaries discussed above. The 
fractional effect of such scattering on the magnon energy 
Eg is on the order of {SL/L)Eq, where 5L is the width of 
the boundary. If we take SL to be one lattice constant, 
the corresponding change in magnon energy is on the or- 
der of 1 /teV at the zone center and 7 fieV at the zone 
boundary. 
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A mechanism that could potentially produce a 
linewidth that increases monotonically along the a-axis 
direction is 4-magnon scattering in which the electronic 
magnons collide with nuclear spin- waves; however, this 
contribution is likely to be more than an order of mag- 
nitude too small[5D]. Other mechanisms which may con- 
tribute, all of which are estimated to be at least an order 
of magnitude too small to explain the data, were sum- 
marized earlier [T^. 

The total calculated linewidth F = r,„ag + agrees 
well with the data (Figs. 2-4). In Rb2MnF4, the 
agreement worsens with increasing T at intermediate q 
(0.025 - 0.15 r.l.u.); in MnFz, for all q, and the dis- 
crepancy is largest at intermediate q (0.2 - 0.4 r.l.u.). 
Additional scattering channels may become important 
at higher temperatures. One candidate is 6-magnon 
(three-in/three-out) scattering, which was posited to pro- 
vide the dominant contribution to the AFMR linewidth 
in high magnetic fields at temperatures above Tjv/2 in 
Rb2MnF4 [21] and MnFj [22]. 

Selected linewidth data from which the fitted curves to 
the 3 K data (Fig. 3) have been subtracted are compared 
in the main panel of Fig. 4b. The exponents from (phe- 
nomenological) power-law fits for these and additional 
values of g/ qzB are shown in the inset. The g-dependence 
of the exponents is non-monotonic; the values for the 
two compounds agree within error at low g, and diverge 
near the respective zone boundaries. The exponent for 
Rb2MnF4 there is 2.6 ± 0.1, which is not too different 
from Kopietz' analytical result of 3 for the isotropic case 
(derived from approximations made to simplify the gen- 
eral expression of TH at the zone boundary) |8j , though 
the region of validity of his approximations lies lower in 
temperature than our data (see Fig. Id). 

Broadening resulting from longitudinal spin 
fluctuations [6] was earlier reported to describe the 
magnon linewidths in MnF2[23]. At the temperatures 
considered [B], the longitudinal fluctuations are in fact 
composed to a large extent of mobile spin- waves, and 
thus collisions with these fluctuations are equivalent to 
the magnon-magnon scattering events calculated in more 
detail here. The self-consistent numerical evaluations 
of scattering within the spin-wave gas performed in the 
current work thus replace both prior comparisons with 
the theory in 3D. 

Comparison of the hydrodynamic prediction with our 
data can be done where sufficient data for a power-law 
fit falls within the low-g, low-T domain of validity of 
the hydrodynamic theory [T]. In our measurements, this 
applies only to MnF2 at 25 K (the upper limit of the 
range of validity), where we obtain an exponent of 2.2 ± 
0.1. Over the range from 15 - 40 K, this exponent ranges 
from 1.9 - 2.2. Hydrodynamics is expected to break down 
in 2D [S]; in Rb2MnF4, the exponent ranges from 1.1 - 
1.3 from 16 - 23 K (however, these temperatures lie above 
the domain of validity). 
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FIG. 4. (Color online) (a) Linewidth data vs. T in Rb2MnF4 
(symbols), measured by NRSE-TAS at wavevectors along the 
{qqQ) direction. The curves are calculated as described in Fig. 
2. (b) Temperature-dependent component of the magnon 
linewidths in Rb2MnF4 (open symbols) and MnF2 (filled sym- 
bols). The fitted contribution from 2- magnon scattering from 
crystal blocks at 3 K (dashed curves from Fig. 3) has been 
subtracted oS. The lines (solid for Rb2MnF4 and dashed for 
MnF2) are power-law fits to the data. The data and curves 
are color-coded by q/qzB- Inset: the exponent corresponding 
to such power-law fits in Rb2MnF4 (red symbols) and MnF2 
(blue symbols). 



In summary, comprehensive comparison of experimen- 
tal data and numerical results suggests that the processes 
that determine the magnon lifetime in 2D and 3D anti- 
ferromagnets are two-in/two-out magnon collisions and 
boundary scattering. The difference between 2D and 3D 
systems is thus quantitative, not qualitative: it origi- 
nates from the difference between near-cylindrical and 
ellipsoidal (near-spherical) surfaces of constant magnon 
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energy with respect to the conservation of energy and 
momentum in 4-magnon collisions. 

The seminal theories of 4-magnon scattering consid- 
ered here could not previously be compared realistically 
and quantitatively with experimental data. This limi- 
tation is relieved here through proper inclusion of the 
Umklapp scattering and the use of numerical evaluations, 
permitting the first complete experimental and theoret- 
ical description of magnon linewidths over broad ranges 
of wavevector and temperature in two and three dimen- 
sions. Comprehensive, momentum-resolved linewidth 
data that are described accurately by theory enable the 
first calculation of the magnon-mediated thermal conduc- 
tivity without adjustable parameters I24j . 
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